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We show that conditional measurements on atoms follow- 
ing their interaction with a resonant cavity field mode can 
be used to effectively counter the decoherence of Fock-state 
superpositions due to cavity leakage. 

PACS numbers: 03.65.Bz, 89.70.+C, 42.50.Dv, 32.80.-t 



I. INTRODUCTION 

Decoherence of non-classical states of a quantum sys- 
tem via coupling to a reservoir is of fundamental interest, 
as it constitutes the mechanism that yields the classical 
limit of the system dynamics jlj . Recently it has also be- 
come a topic of great applied importance, because it de- 
termines the feasibility of quantum information storage, 
encoding (encrypting) and computing In numerous 
current theoretical proposals, the irreversibility of deco- 
herence processes in quantum computing is combatted by 
two generic means. One is the filtering out of the ensem- 
ble portion which has not decohered, i.e., has remained 
intact. This approach has been suggested for two-mode 
fields ||, but not for single-mode cavity fields. The other 
means is encoding the state (qubit) by means of several 
ancillas, decoding the result after a certain time, checking 
the ancillas for error syndromes and correcting them . 
Although the latter approach is in principle applicable 
to arbitrary errors, only extremely small error probabil- 
ities (per qubit or gate, per time step) can afford fault 
tolerant quantum computation Q . Instead of the "high 
level" unitary transformation approach to error correc- 
tion in quantum computing — which involves substantial 
overhead in qubits and gates — the countering of deco- 
herence of stored quantum information, e.g., in between 
computation steps, may be achieved by a "low level" ap- 
proach: applying simple physical manipulations to the 
quantum storage device, which take advantage of its spe- 
cific physical realization. Such approach has been advo- 
cated recently ||, relying on continuous monitoring of 
the dissipation channel for quantum jumps, with perfect 
photodetection efficiency, and on instantaneous feedback 
for the inversion of their effect. 

The nature of quantum computing requires that deco- 
herence be corrected without knowing which state is in 
error during the computation. There is, however, a sim- 
pler but still important problem: how to protect from 
decoherence the input states, prior to the onset of com- 
putation. Here we suggest a non-unitary approach to 



counter decoherence, which can be used to safely store 
quantum field states in dissipative cavities, in order to 
subsequently use them as apriori known input in informa- 
tion processing or in signal transmission. The basic idea 
is to restore the decohered field state by entangling it with 
an atom, and then projecting the entangled state onto a 
superposition of atomic eigenstates, whose phase and am- 
plitudes are specifically tailored for the field state we wish 
to recreate. Such projection amounts to post-selection of 
the appropriate atomic state, i.e., to a conditional mea- 
surement (CM) @J^]. The specific scheme we put forward 
is based on modification of our optimized CM strategy 
Q for cavity-mode state preparation by resonant inter- 
action with atoms, in the Jaynes-Cummings (JC) model, 
followed by projection onto selected atomic states. In the 
present problem we set the initial (unspoilt) superposi- 
tion of zero-photon up to A^-photon states as our tar- 
get state, and work in Liouville space instead of Hilbert 
space, so as to account for the state decoherence. The 
results demonstrate that a few highly-probable CMs, in 
this simple model, can drastically reduce even a large er- 
ror. One of our objectives is to find the optimal tradeoff 
between the CM probability and the error size, which 
grows in the course of dissipation. 

The ability to approximately restore any mixture to 
any pure state (in our A^+l-dimensional Hilbert space) is 
the advantage of our post-selection CM approach, com- 
pared to the non-selective measurement (tracing) ap- 
proach: Mixed states can only evolve into the special 
"cotangent" and "tangent" pure states by a large num- 
ber of JC interactions with atoms initially prepared in su- 
perposition states (under the atomic excitation-trapping 
condition) followed by tracing over the atomic states jio) . 



II. DECOHERENCE MINIMIZATION BY 
CONDITIONAL MEASUREMENTS 

We consider a single-mode cavity in which the quan- 
tized electromagnetic field is initially prepared in a finite 
superposition of Fock states, 
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To model the effect of dissipation we assume the cavity 
field to be coupled to a zero-temperature heat bath. The 
master equation describing such coupling, in the interac- 
tion picture, is 
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where p F = p F (t) is the density matrix of the cavity field, 
a and a) are the annihilation and creation operators of 
the field, and 7 is the damping constant of the cavity. 

The solution of Eq. (Q) after dissipation over time i > 
[fill can be shown to have the form 
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written here in Fock basis, p n , m (t) = (n\p F (t)\m). 

In order to recover the original state of the field we pro- 
pose to apply an optimized CM (or a sequence thereof) 
to the cavity as follows: Using a classical field we prepare 
a two-level atom in a chosen superposition |9|Jl2]] 
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of its ground \g) and excited |e) states, and let it in- 
teract with the field for a time r by sending it through 
the cavity with controlled speed. The field-atom inter- 
action is adequately described by the resonant Jaynes- 
Cummings (JC) model 0|. We assume the field-atom 
interaction time r to be much shorter than the cavity life- 
time, 7T -C 1, so that we may neglect dissipation during 
each CM. Upon exiting the cavity the atom is condition- 
ally measured, using a second classical field, to be in a 
state 
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differing in general from the initial atomic state \4>^)- 
This means that we post-select, using the same setup 
as in ref. ||, the atomic superposition state (||) which is 
correlated to a cavity field state that is as close as possible 
to the original state (A3). 

The effect the applied CM has on the cavity field is 
then calculated as follows: Initially, at the time the atom 
enters the cavity, the density matrix of the field-atom 
system is 
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It then evolves unitarily by the JC interaction of duration 
t into 



p FA (t + r) = U(r)p FA (t)UHr) 
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where U(r) is the interaction picture evolution operator 



U{T)\n)\e)=C n \n)\e)-iS n \n + l)\g) 
U{r)\n)\g) = C„_i|n)|<?) - iS n ^\n- 1) 



(8) 



with C n = cos (AtV n + l) and S n — sin [Xr^/n + l) , A 
being the field-atom coupling constant (known as the vac- 
uum Rabi frequency). Finally, the conditional measure- 
ment of the atom in the state \<j>^) results in a density 
matrix of the field given by 
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is the success probability of the CM. The explicit expres- 
sions for p F (t + t) and P are given in the Appendix for 
an initial superposition of |0) and |1) states. 

To nearly recover the original state of the field, we use 
the dependence of p F (i+T) on the initial and final atomic 
states and the field-atom interaction time, choosing op- 
timal parameters aW, j3^ % \ fjw and r such that 
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holds (see Appendix for an explicit form of this con- 
dition), along with high CM success probability (|ic|). 
These optimal CM parameters are found by minimizing 
the cost function || 



G 



d(p F {i+T) lPF (0)) 
pr 



(12) 



where d is a distance function between two density ma- 
trices, defined as 
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P is the CM success probability (|l0|), and the adjustable 
exponent r > determines the relative importance of the 
two factors in G. If this CM does not bring us as close to 
the original state as our experimental accuracy permits, 
we can repeat the process over and over again, as long as 
the distance to the original state keeps decreasing, while 
the CM success probability remains high. The atomic 
states (|J) and (^|) are determined by the minimization of 
( |l2| ) at each step. Let us note here that the application of 
each CM may introduce widening of the photon-number 
distribution by one photon, and yet the optimized CMs 
are capable of avoiding this widening and, moreo ver, o f 
resto ring the field to its initial pure state. Eqs. ( AlOe - 
AlOi ) in the Appendix exemplify the widening-avoidance 
requirements which are implicit in condition ([ll]) . These 
requirements amount to an effective control of a large 
Fock-state subspace. 



III. EXAMPLES 

We illustrate our approach with two examples below, 
using the Q-function Q Pp (a, a*) = (a\p F \a), \a) being a 
coherent state of complex amplitude a, to visualize the 
error-correction process: 

1) Let us take as the original field state an equal- 
amplitude superposition of our basis states, e.g., 
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whose Q-function is shown in Fig. 1(a). Dissipation 
by yt = 0.3 renders the error matrix p F (t) — p F (0) 
of considerable magnitude, as seen in Fig. 1(b). Af- 
ter the application of one CM (|</>«) = cos(37r/8)|e) + 
sin(37r/8)e l57r / 4 | 5 ), At = 37.95, |<^) = cos(37r/8)|e) + 
sin(37r/8)e lTr / 4 |(7)), optimized to yield high success proba- 
bility (r = 2), the remaining error matrix p F (t+r)—p F (0) 
is roughly 2.5 times smaller than before the correction, 
as seen in Fig. 1(c). The success probability of the CM 
is a high 74%. Subsequent CMs can further reduce the 
distance to 1/6 (one sixth) its original magnitude, with 
62% success probability for the full CM sequence (Fig. 3). 
Stronger error reduction is obtainable at the expense of 
success probability: the application of 4 CMs optimized 
for r = 1 (respectively r = 0) yields an error reduction 
factor of 11 (respectively 28) with sequence probability 
of 33% (respectively 16%). 

2) If the original field state is a strongly unequal su- 
perposition of the basis states, such as 
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(Fig. 2(a)), the error matrix after dissipation by ji = 0.3 
is again significant (Fig. 2(b)). Successive application of 
4 CMs, optimized for r = 2, reduces this error by a factor 
of 30 (Fig. 2(c)), which means that the recovered state is 
practically indistinguishable from the original state. The 
success probability of the total CM sequence, 50%, is 
markedly high (Fig. 3). If we ignore success probability 
in Eq. (12) (r = 0) we obtain a higher error-reduction 
factor of 75, with sequence probability of 28%. 

In Fig. 3 we plot the distance dx = d(p K , p F (0j) (Eq. 



(|13|)) between the recovered state and the original state 

and the CM sequence probability P seq ,K — Yii=i Pii with 
Pi given by (|l0|), as a function of the number of CMs 
performed. It shows that the first CMs achieve a strong 
reduction of such a distance, whereas after a few succes- 
sive CMs saturation sets on, in terms of both distance 
and success probability. 

It is interesting to compare the success probability 
in our approach with the theoretical probability to find 
the original state in the dissipation-spoilt state, namely, 
Tr F [p F (0)p F (t)], which we call the filtering probability. 
In Table I we list the success probability of a sequence 
of 4 CMs (optimized for r = 2), P seq ,K=i> and the cor- 
responding filtering probability for various values of the 
dissipation parameter ji, taking as the original state the 
state © of example 2. The probability Tr F [p F (0)p* = 4 ] 
of finding the original state in the recovered state is 0.99 
or higher for all entries. 



IV. DISCUSSION 

In conclusion, we have demonstrated here the effective- 
ness of simple JC-dynamics CMs as a means of revers- 
ing the effect of dissipation on coherent superpositions of 
Fock-states of a cavity field: the application of a small 



number of optimized CMs recovers the original state of 
the field with high success probability, which is com- 
parable or even surpasses the filtering probability. The 
simplest tactics may employ a single highly-probable trial 
to achieve nearly-complete error correction. As noted 
above, although we have only five control parameters at 
our disposal for each CM, our optimization procedure is 
able to effectively control the amplitudes in a large Fock- 
state subspace. 

Among the experimental imperfections that can de- 
grade the effectiveness of any CM approach j^,[l4|,|l5| , re- 
alistic atomic velocity fluctuations (of 1%) and cavity- 
temperature effects (below 1°K) are relatively unimpor- 
tant, and especially so in the present scheme which makes 
use of a single or few CMs so that the effect of ex- 
perimental imperfections is linear in the input errors. 
Only atomic detection efficiency is an experimental chal- 
lenge |lj] . Although the detection efficiency is currently 
low, it is expected to rise considerably in the coming fu- 
ture. 

Extensions of this approach to field-atom interaction 
Hamiltonians with more controllable degrees of free- 
dom fll4| can make a single trial within this correction 
procedure effective for highly complicated states, encod- 
ing many qubits of information. Nevertheless, even in 
its present form the suggested approach has undoubted 
merits: (a) it can yield higher success probabilities than 
the filtering approach; (b) it is not limited to small er- 
rors as "high level" unitary-transformation approaches 
are; (c) it corrects errors after their occurrence, with no 
reliance on ideal continuous monitoring of the dissipation 
channel and on instantaneous feedback; and (d) it is re- 
alistic in that it can counter combined phase-amplitude 
errors which arise in cavity dissipation, and is of general 
applicability — not restricted to specific models of dissi- 
pation. 
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APPENDIX: 

The reduced density matrix of the field resulting from 
its interaction with the atom followed by the conditional 
measurement on the latter can be found using the for- 
mula (Eqs. (H0)) 
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where the normalization constant P is the success prob- 
ability of the conditional measurement and is given by 



P = Tr„ Tr 



(A2) 



In the simple case where the initial field state is a super- 
position of the vacuum and one-photon states 

^(0))=co|0) + Cl |l), MO) = |#))><V>(0)| , (A3) 

the density matrix resulting from dissipation over time t 
is 
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The explicit expressions for p F (t + r) and the success 
probability P are then 

p F (i+r) = P~ l { [\a^\ 2 A + a W0V>M + ^a ({ >C 
+ |/3 (f) | 2 0] |0) (0| + \\a^\ 2 B + a (f) /? (f) *A + (3 {i) a^*D 
+ \/3^\ 2 k] |0)(1| + \\a^\ 2 E + a (f) /? (f) *i? + f3^a {i) *H 
+ |/3( f )| 2 rl |1)(0| + [|a (f) | 2 F + a®p®*S + p {i) a^*I 



|1)(1| + \a ( - f) [3 ( - l) *W +\(3 (l) \ 2 Y] |2)(0| 



(f)* G+ |^(f)|2gj | > <2| + [a (f) /3 (f) *A 
|2)<l| + [/3< f V f )*L+|/3»|V] |1)(2| 
|2)(2|}, 



and 



P = \a^\ 2 (A + F) + a^(3 < - l) *(M + S) 

+ /? (f) a (f) *(C + 1) + |/3 (f) \ 2 {0 + U + Z) 

The coefficients A,B,... here are given by 

A = poo(t)\a^\ 2 C 2 +ipoi(t)a^p^*C So 

+ p n (i)\f3 {i) \ 2 So ~ zpio SqCo , 
B = poi^laUfCod - ipuffiaP'pVdSo , 
C = poo(i)««/3«*C - z P io(t)|/3 (i) | 2 So , 
D = ip o(t)\a^\ 2 C a So - t Pll (t)\l3 (i) \ 2 S Q C 
+ P01 (t)a«/?«*C 2 + p 10 (i> (i) */? (i) S 2 , 
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F = Pll {t)\a^\ 2 Cl , (A8e) 

G = pnffiaV'flVSoSi + i Pm {t)\a^ | 2 C ^i , (A8f) 

H = p lo {i)a®0®*C 1 , (A8g) 

/ = pii(t)aW/3«*aiC + ipio(t)|a (i) | a <7iSb , (A8h) 

X = ip o(t)a (i) */3 (i) ^o + Poi(t)|/3 (i) | 2 C , (A8i) 

i = i Pll (t)|a«| 2 C 1 5 1 , (A8j) 

O = poo(i)|/3 (i) | 2 , (A8k) 

Q = ip i(t)a (i) */3 (i) ^i , (A81) 
C/ = p Q0 (t)\a^\ 2 S 2 - ipoi®aF>0®*SoC o 

+ Pll (t)\P^\ 2 C 2 +ip w (F) a < - i >^C So , (A8m) 

K = ipufta^pVCoS! +p i(t)\a^\ 2 S S 1 , (A8n) 

Z = m (t)|a«| 2 5 2 , (A8o) 

with the following relations holding between them 

A = A* , B = E* , C = M* , D = R* (A9a) 

F = F* , G = W* , H = N* , I = S* (A9b) 

L = A* , = 0*, K = T* , Q = Y* (A9c) 

U = U* , V = J* , Z = Z* . (A9d) 

(The coefficients A and A in this appendix bear no re- 
lation to A and A mentioned in the main text). 

The explicit form of condition (|ll] ) for recovering the 
original field state is given by the following list of approx- 
imation relations: 
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FIG. 1. Q-function description of (a) original field, with 
M0) = |V(0)>(lK0)|, |V>(0)) = (|0) +e^ 3 |l))/V2; (b) er- 
ror after dissipation, p F (t) — p F (0) ; (c) reduced error after 1 
optimized CM (minimizing Eq. (J12|)) , p F (i + r) — p F (0). 



5 




FIG. 2. Q-function description of (a) original field, 

p F (0) = |V(o))(V(o)|, |V(o)> = 10^10) + e ^Vi-io- 2 |i) ; 

(b) error after dissipation, p F (t) — p F (0); (c) error after 4 
consecutive optimized CMs, each minimizing Eq. jl^). 



K 

FIG. 3. Distance dn and probability P S eq,K as functions of 
the number K of CMs applied to the cavity. Solid curves - 
example 1; dashed curves - example 2. 



TABLE I. Comparison of CM success probability with fil- 
tering probability. 



7* 


0.3 


0.4 


0.5 


1.0 


Pseq,K—A 


0.50 


0.35 


0.40 


0.63 


Tr F [p F (0)p F (t)] 


0.56 


0.46 


0.38 


0.15 
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